Strongly anisotropic roughness in surfaces driven by an oblique particle flux 
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Using field theoretic renormalization, an MBE-type growth process with an obliquely incident 
influx of atoms is examined. The projection of the beam on the substrate plane selects a "parallel" 
direction, with rotational invariance restricted to the transverse directions. Depending on the behav- 
ior of an effective anisotropic surface tension, a line of second order transitions is identified, as well 
as a line of potentially first order transitions, joined by a multicritical point. Near the second order 
transitions and the multicritical point, the surface roughness is strongly anisotropic. Four different 
roughness exponents are introduced and computed, describing the surface in different directions, in 
real or momentum space. The results presented challenge an earlier study of the multicritical point. 
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I. INTRODUCTION 

The fabrication of numerous nanoscale heterostruc- 
tures requires the controlled deposition of material onto 
a substrate. A variety of deposition processes are used, 
depending on desired surface structure and device per- 
formance. Molecular beam epitaxy (MBE), involving di- 
rected beams of incident atoms, is particularly suitable if 
lower growth temperatures and precise in situ control and 
characterization are desired^. It is an important goal of 
both theoretical and experimental investigations to gain 
an understanding of the resulting surface morphology, in 
terms of its spatial and dynamic height-height correla- 
tions, or more specifically, its roughness. 

Beyond the obvious implications for nanoscale devices, 
surface growth problems also constitute an important 
class of generic nonequilibrium phenomenal. Particles 
are deposited on the surface and may diffuse around 
on it. If deposition occurs from a vapor, desorption or 
bulk defect formation tend to be important processes; 
in contrast, both mechanisms can often be neglected in 
MBE (see for example^). After an initial transient, a 
steady state is established which is characterized by time- 
independent macroscopic properties, provided a suitable 
reference frame is chosen. Generically, detailed balance is 
broken by the incident particle flujs^, so that this steady 
state cannot be described by a Boltzmann distribution; 
instead, its statistical properties have to be determined 
directly from its dynamical evolution. If one is primar- 
ily interested in universal, large scale, long time charac- 
teristics, the dynamical evolution can often be cast as 
a Langevin equation which can be analyzed using tech- 
niques from renormalized field theory. 

Here, we extend a model^ due to Marsili et al. to de- 
scribe MBE-type or ballistic deposition processes with 



oblique particle incidence. Focusing on large scale prop- 
erties such as surface roughness, we exploit a coarse- 
grained (continuum) approach. Adopting an idealized 
descriptionSiSiS, particle desorption and bulk defect for- 
mation will be neglected so that all (deterministic) sur- 
face relaxation processes are mass-conserving, i.e., can 
be written as the gradients of equilibrium and non- 
equilibrium currents. Shot noise in the deposition process 
requires the addition of a stochastic term to the growth 
equation. Since the particle beam selects a preferred 
("parallel") direction in the substrate plane, the resulting 
Langevin equation is necessarily anisotropic. The inter- 
play of interatomic interactions and kinetic effects, such 
as Schwoebel barriers, generates an anisotropic effective 
surface tension which can become very small or even van- 
ish. Due to the anisotropy, this leads to four different 
regimes with potentially scale-invariant behaviors. We 
analyze these four regimes, identify the scale-invariant 
ones, and compute the associated anisotropic roughness 
exponents. 

Models with oblique particle incidence have been in- 
vestigated previously. Focusing on vapor-deposited thin 
films, Meakin and KrugiSiiiii^ considered the ballistic de- 
position of particles under near-grazing incidence. Under 
these conditions, columnar patterns form which shield 
parts of the growing surface from incoming particles. The 
large scale properties of these structures can be charac- 
terized in terms of anisotropic scaling exponents, differ- 
entiating parallel and transverse directionsiSii^. 

Following Marsili et ali^, our model differs from Krug 
and Meakin's approach in two important respects. First, 
surface overhangs and shadowing effects are neglected 
so that our results are restricted to near-normal in- 
cidence. Second, our model is designed for "ideal 
MBE" -type growth so that it falls outside the Kardar- 



Parisi-Zhang universality classic for mass non-conserving 
growth. However, extending the work of Marsih et aJ^, 
we include a possibly anisotropic effective surface tension 
and investigate its effects systematically (see further com- 
ments in the following section). Due to the anisotropy, 
this contribution actually consists of two terms, one con- 
trolling relaxation parallel and the other transverse to 
the beam direction, with coupling constants T|| and r^, 
respectively. Depending on which of these two couplings, 
T|| or T_L, vanishes first, ripple-like surface structures are 
expected, aligned transverse to the soft direction. The 
roughness properties near these two instabilities, charac- 
terized respectively by r_L = while T|| > 0, and t± > 
while T|| =0, are discussed in this paper for the first time. 

The original theory of Marsili et al^ is recovered only if 
both couplings, t\\ and t_l, vanish simultaneously. Since 
the latter requires the careful tuning of two parameters, 
it is much less likely to be experimentally relevant than 
either of the other two instabilities in which only a single 
parameter must be adjusted. Referring all technical de- 
tails to a separate publication^^, we will point out very 
briefly that the fixed point and the roughness exponents 
reported in Marsili et a& must be significantly revised. 

An important aspect discussed in this article is how 
to translate surface data into roughness exponents, for 
inherently anisotropic surface models such as the one an- 
alyzed here. If we generalize the standard definitions 
familiar from isotropic problems, we arrive at four dif- 
ferent roughness exponentsii^ Two of these characterize 
real-space scans along and transverse to the beam di- 
rection, and the remaining two are needed to describe 
scattering (i.e. momentum space) data with parallel or 
transverse momentum transfer. All four are related by 
simple scaling laws, but possess distinct numerical val- 
ues. When analyzing experimental data, it is therefore 
essential to be aware of these subtleties. 

The paper is organized as follows. We first present the 
underlying Langevin equation for a single-valued height 
field and briefly review the physical origin of its con- 
stituents. We then present a careful definition of the 
roughness exponents, based on height-height correlation 
functions and their structure factors. Turning to the 
renormalization group (RG) analysis, we first discuss a 
simple scaling symmetry of our model which allows us to 
identify a set of effective coupling constants. The invari- 
ance of our model with respect to tilts of the surface is 
much more powerful. Since this symmetry is continuous, 
it gives rise to a Ward identity which relates different ver- 
tex functions. This simplifies the renormalization proce- 
dure considerably. We then present our main results for 
the scaling properties of correlation and response func- 
tions for the four different cases: (o) both rii and tj_ are 
positive; (i) rii remains positive while t± vanishes; (ii) tu 
vanishes while t± remains positive; and finally, (iii) both 
T|| and T_L vanish simultaneously. Roughness and dy- 
namic critical exponents are derived. We conclude with 
a short summary and a discussion of the experimental 
evidence. 




Figure 1 : Cartoon of the arrangement of surface and incom- 
ing flux. The field h,(x, t) is the height of the interface over 
the substrate. Its coordinate system can always be chosen 
such that a; II is parallel to the projection of the flux on the 
substrate. The remaining component of the latter is the (neg- 
ative) flux Jz- 



II. THE MODEL 

We focus on long time, large distance phenomena of 
the growing surface. Under suitable condition^, sur- 
face overhangs and shadowing effects may be neglected, 
so that the surface can be described by a single-valued 
height field, h{r, t), where r denotes a d-dimensional vec- 
tor in a reference (substrate) plane, the z-axis is normal 
to that plane, and t denotes time, see Fig. ^ The no- 
overhang assumption can be justified a posteriori if the 
calculated interface roughness exponent are found to be 
less than unity. The time evolution of the interface is 
described by a Langevin equation of the form 



dth = G[h] + r/ 



(1) 



Here, r/ denotes the effects of shot noise, and G[h] models 
the deterministic part of the surface evolution, assumed 
to be mass-conserving in a suitable coordinate system so 
that G[h] can be written as a divergence, G[h] = V • F[h]. 
One contribution to G[h] is due to the incident flux; the 
other contribution arises from surface diffusion. All of 
these contributions can be derived using the principle of 
re-parametrization invarianco^i and arc discussed in the 
following. 

As shown in Fig. ^ the incident particle current has 
a normal component Jz and a component parallel to the 
substrate plane, J||. In the following, the co-ordinate 
system is rotated such that one of the axes, labelled a;|| , is 
aligned with Jn . The particles themselves arc of finite size 
with radius r,,, which implies that the flux responsible for 
growth at some point on the surface is to be measured at 
a distance r normal to the surface, see Fig.|3 This effect 
has been discussed in detail in the literature, seoiiiSiiS. 
Neglecting higher order terms in the spirit of a gradient 
expansion, to leading order, this effect gives rise to a 
deterministic term of the form 

Gdrivc[h] = -J, + J|| • V/i + roJ^V'^h - ra (J|| • Vh) V^/i 

(2) 
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Figure 2: Schematic representation of the growth inhibition 
and amplification by either finite adatom radii (ro) or steering. 
In both cases, the fiux contributing to growth at some point 
on the surface is to be measured a distance To normal to the 
surface. The left box shows the reduced effective surface in 
valleys, the right one the corresponding effect on peaks—*—. 



Most remarkably, so called "steering" leads to the same 
terms in leading order—. Steering implies that deposited 
atoms are deflected towards the surface normal as soon 
as they reach a certain distance above it, due to an at- 
tractive force exerted by the particles in the deposit. 
Returning to Eq. ^^ we note that the first two terms 
in Gdrive[/i] can be removed by a Galilei transformation 
h{r,t) — > /i(r + J||t, i) — J^t. From now on, we always 
work in this co-moving frame. 

In addition to being driven by the incident flux, 
Eq. |(2Jl, the surface also relaxes via diffusion along the 
surface, leading to a quartic term of the form^ 



G 



relaxation 



= -^V^/i 



(3) 



In the following, terms of the form rV^h play a partic- 
ularly important role, since they determine which of the 
various critical regimes can be accessed. Eq. Q already 
contains such a term, induced by Jz, but other contri- 
butions of this type are possible, e.g., a negative term 
from a step edge (Schwoebel) barrieriSSiS^ or a positive 
one due to a surface tension. Moreover, even if such a 
term were initially absent, it would actually be gener- 
ated under renormalization group transformations and is 
therefore intrinsically present. In contrast to Marsili et 
al^, we include it from the very beginning. 

It is essential to note that the nonlinear term in Eq. ((^J, 



(J|| • \/h) 



V^/iEE A(5||/i) V^/i 



introduces an anisotropy into the system which breaks 
the full rotational symmetry within the d-dimensional 
space of the substrate. As a consequence, there is no rea- 
son to expect isotropic coupling constants (such as /i or t) 
for the linear contributions. Instead, any coarse-graining 
of the microscopic (atomic level) theory is expected to 
give rise to different couplings (such as tii and r^) in 
the parallel and perpendicular subspaces, and this is in- 
deed confirmed by the renormalization group. If these 
anisotropics are incorporated into the theory, preserv- 
ing only rotational invariance in the {d — l)-dimensional 
transverse subspace, the Langevin equation ^ takes the 
form 



- pi^{V\fh + {d\\h) {\\dlh + \x_V\h) + ry 

= -a||j|| - V Ji + 7/ (4) 

where we have introduced an explicit time scale for con- 
venience. The surface currents are given by 

-Ai(Vi/i)(9||/i). (5) 

The scalar differential operator dn operates only along x\\ 
(see Fig.^, while the vector Vj_ operates in the (d — 1) 
dimensional subspace perpendicular to x\\. We also note 
that the nonlinearity now splits into two distinct terms, 
with couplings Ay and A^, respectively. Finally, the con- 
served nature of the deterministic surface evolution is 
displayed explicitly here. 

The randomness of particle aggregation on the surface 
is captured by the non-conserved white noise ri[r, t) with 
zero average and second moment 



(77(r,t),?(r',0>=27-i<5(r-r')5(t-t') 



(6) 



Eq. I0J forms the basis for the following analysis. Its 
properties are controlled by the dominant terms in the 
gradient expansion. To ensure the stability of the linear 
theory, we require ^^\\,^J.A_ > 0, and ^x > — (m||M-l)^ ■ 
The two couplings tk and t±_ play the role of critical 
control parameters. If both are positive, the nonlinear- 
ities become irrelevant, and the problem reduces to the 
well-known Edwards- Wilkinson equation^. In contrast, 
if one or both of them vanish, the long time, long dis- 
tance properties of the theory change dramatically: The 
surface undergoes an instability and forms characteristic 
spatial patterns. If just one of the couplings goes soft, 
these patterns take the form of ripples (similar to corru- 
gated roofing) transverse to the soft direction. If both 
couplings become negative, the surface develops mounds 
or "wedding cakes"—. Focusing only on the onset of these 
instabilities, four different cases emerge whose properties 
are discussed in the following: (o) the "disordered" phase, 
corresponding to the linear theory with rii > 0, r^ > 0; 
(i) a line of continuous transitions t\\ > 0, r^ ^ 0; (ii) 
a line of possibly first order transitions r|| ^ 0, tj^ > 0; 
and (iii) the multicritical (critical end-) point rii — > 0, 
Tj^ -^ 0. Before we turn to any technicalities, however, 
we first discuss an important physical issue, namely, the 
definition of appropriate roughness exponents. 



III. ANISOTROPIC ROUGHNESS EXPONENTS 

The roughness of the surface, and the associated rough- 
ness exponents, if they exist, are easily measured exper- 



imcntally. They can be determined from real-space im- 
ages of the surface or from scattering data in momentum 
space. Within our theoretical framework, roughness ex- 
ponents can be extracted from the height-height correla- 
tion function, 



C(r-r',t-t')= (/i(r,t)/i(r',t')) 



(7) 



Since we focus on the steady state in the absence of spa- 
tial boundaries, we assume translational invariance in 
space and time. The spatial Fourier transform of C is 
the dynamic structure factor, 



C{ci,t)^Jd''rC{r,t)e^^' 



In the absence of anisotropics, the asymptotic scaling 
behavior of Eq. jT)) can be written in the form 



C{r,t)^\r\'^c{t/\rn 



(8) 



where x denotes the roughness exponent and z the dy- 
namic exponent of the surface while c is a universal scal- 
ing function. In Fourier space, the behavior of C(r, t) 
translates into 



C{q,t) 



-{d+2x) 



c{t\q[ 



(9) 



In the presence of strong anisotropy, where the scaling 
of the correlation function depends on the direction, the 
situation becomes more complex. Anticipating some of 
the following results, a key finding of the present arti- 
cle is the existence of four different roughness exponents, 
characterizing the surface along the parallel or transverse 
directions, in real or in momentum space. While they 
are directly related by scaling laws, it is essential to re- 
alize that they take different numerical values. In the 
interpretation of actual experimental data, it is therefore 
important to identify the appropriate member of this set 
of four exponents in order to compare with theoretical 
predictions. 

In our scaling analysis below, we adopt the conven- 
tional exponent definitions from critical dynamics. The 
exponent v controls the divergence of the correlation 
length, while 77 denotes the anomalous dimension of the 
height field which appears in all correlation functions. 
Due to the presence of anisotropy, an additional expo- 
nent, i.e., the strong anisotropy exponent A, is required 
to reflect the different scaling of distances or wave vec- 
tors in different directions'^. If I denotes an arbitrary 
transverse momentum scale, so that |q_L| ex /, we intro- 
duce A via q\\ ex /^+^. With this definition, one finds 
that, asymptotically, the structure factor is a generalized 
homogeneous function of its variables 

Ciqpqr.t) - r4+"C(<Z||//l+^,qx//,Pt) (10) 

and in real space, the two-point function takes the form 

C{x\\,Yx^;t) ^ V^+^-'^+'^Cif^+^X\\jY^,lH) (11) 



In analogy to Eq. © two roughness exponents are defined 
in real space, Xi- and X||j via 



C(0,rx;t) 

C(x||,0;t) 



J^^-c^(t/|r^r) 



|X||| "C||(Va: 



2/(1+A)n 

II > 



(12) 



Of course, this is only meaningful if the two scaling func- 
tions Ci_ and c approach finite and non-zero constants 
when their arguments vanish. Under this assumption, 
the two exponents 

X± = ^[4-(d+A)-,7] 

XII = i(l + A)-i[4-(rf + A)-r7] (13) 

are read off immediately. In order to define the corre- 
sponding exponents in momentum space, xj_ and X|| , we 
focus on two structure factors which are easily accessible 
experimentally, especially if we set t = Q: 



C(0,qi,i) =|qx 



~{d+2xi_) 



c_L(i|q±| 



C(g||,0,i) =q. 



_ -(d+2x||)- 



i9| 



II 



(14) 



Again, provided the scaling functions cx and c\\ are non- 
singular and non-zero in the limit of vanishing argument, 
one reads off 



1 



X± 
XII 



^ 



2^1 + A J 



(15) 



The key observation is that x±_ = X-\- and xii = x\\ only 
if the anisotropy exponent A vanishes. Thus, contrary to 
Eqs. ((SJ and ^ which are equivalent definitions of the 
roughness exponent in isotropic systems, the correspond- 
ing definitions for anisotropic systems will typically give 
rise to different exponents. 

In the following, we explicitly compute the scaling ex- 
ponents in the previous expressions, and also confirm 
the underlying scaling form. To unify the discussion, 
we first recast the Langevin equation as a dynamic field 
theory, collect the elements of perturbation theory, and 
then identify the upper critical dimensions and marginal 
nonlinearities for the four cases. Our final goal is a sys- 
tematic derivation of the scaling properties of correlation 
and response functions. 



IV. RENORMALIZATION GROUP ANALYSIS 

A. Povi^er counting and mean-field exponents. 

In this section, we assemble the basic components of 
the field-theoretic analysis, leaving technical details toi^. 
The formalism becomes most elegant if we introduce a 



response field h{r, t) and recast the Langevin equation 
Q as a dynamic functional J'[h,h], following standard 
methodsS^iSSiSi: 

J[h,h]^j f d'^xdtihU-^dth + d^^^ + VJ_l] -h^Y 

(16) 
This has the advantage that both correlation and re- 
sponse functions can be computed as appropriate func- 
tional averages, with statistical weight exp(— J7). The 
analysis can be simplified considerably if we exploit 
the symmetries of Jj[h^h], using the explicit forms of 
the currents, Eq. (0. First, the symmetry /i(r,t) — > 
h{r,t) + a implies invariance under a coordinate shift 
in the z-direction. Second, the theory is invariant un- 
der tilts of the surface by an infinitesimal "angle" b, 
i.e., h{r,t) -^ h{r,t) -I- b • r provided the tilt is accom- 
panied by a transformation of the couplings, namely. 
Til -^ Til — 6j|A|| and tj_ ^ tj_ — bnX±. Third, parti- 
cle conservation on the surface leads to invariance un- 
der the symmetry transformation h{r,t) — > h(r,t) + c, 
h{r,t) -^ h{r,t) + 2c'yt. Finally, we have a symmetry 
under inversion, namely, h{x\i,rj_,t) — > —h{—x\i,r±,t), 
h{x\\,r±,t) — > —h{—x\\,r±,t). The most important of 
these symmetries is the tilt invariance. Thanks to the 
associated Ward-Takahashi identitjiSSiSS, the renormal- 
izations of rii and t± can be related to those of An and 



Aj^, so that some exponent relations will be valid to all 
orders in perturbation theorjii^. 

Due to the anisotropy, there are two independent 
length scales. If only parallel lengths are rescaled, via 
Xi\ — *■ ax\i, the functional remains invariant provided 
h -^ a^^/'^h, h -^ a^^/'^h and /i|| -^ Qf''^||, T|| — > a^T||, 

^x -^ a^fJ-x while A|| —> a'^^^\\\ and A^ —> a'^/^Ax- 
Likewise, if only transverse lengths are rescaled, via 
r_L —f /3rx, the functional remains invariant provided 
h -^ /3-('^-i)/2/i, h ~> /3-('^-i'/2/i and fij^ ^ /3V±, 



/3^TJ_, ^x 



/32^x while All ^ /3(''-i)/2A|| and 



T± -^ 

A^ -^ /3'^''+'^^/^A_L. As a result, the theory naturally gives 
rise to effective expansion parameters which are invari- 
ant under these rescalings. The precise forms of these 
parameters differ slightly for the four cases and will be 
discussed next. 

In addition to these purely spatial rescalings, we can 
perform a more general dimensional analysis of Eq. (|16|l . 
involving both spatial and temporal degrees of freedom. 
It is standard to express it in terms of an external length 
scale K~^. The key to the scaling of different terms in the 
functional lies in the behavior of the control parameters 
Til and Tj_. Depending on whether they vanish or remain 
finite, the Gaussian part of the dynamic functional is 
dominated by different terms in the (infrared) limit of 
small momenta and frequencies. 

Case o: If both T|| and tj_ are finite and positive, the 
theory turns out to be purely Gaussian. Quartic deriva- 
tives can be neglected in the infrared limit. It is nat- 
ural to scale both parallel and transverse momenta by 
K, via |q_L| oc g|| ex k. As a result, time scales as k~^ 



and the fields have dimensions h{r,t) ex n^'^ ^^/^ and 
h{r,t) ex n'-'^+^y^ so that the nonlinear couplings scale 
as A|| ex A^ ex k^'^''^ and are therefore irrelevant in 
any dimension d > 0. The resulting theory is a sim- 
ple anisotropic generalization of the Edwards- Wilkinson 
equatio] 



.24 



j-^dfh = T\id(\h + Tj_V'_ 



V 



The anisotropics in the quadratic terms affect only 
nonuniversal amplitudes and can be removed by a simple 
rescaling, without losing any information of interest. As 
is well known, the two-point correlation function scales 
as 



C(r,f) 



Irl'-'cit/lrl') 



from which one immediately reads off the (isotropic) 
roughness exponent x = (2 — d,)/2 and the dynamic ex- 
ponent z = 2. Since this case is so familiar (seo^ for a 
detailed discussion), it does not need to be considered 
further. 

Case i: If T|| remains finite and positive but t± van- 
ishes, the two leading (Gaussian) terms in the dynamic 
functional arc ii±h{'V']_)'^h and ^hduh. This suggests 
that parallel and transverse momenta scale differently, 
already at the Gaussian level, namely, |q^| ex k and 
q|| ex K^. If we rewrite the scaling of parallel mo- 
menta as (7|| ex k^"*" , the anisotropic scaling exponent 
A equals unity for the Gaussian theory. Time scales 
as K~^, and t± ex k^ is strongly relevant. fi± can be 
set to 1 by a transverse rescaling with an appropriate 
/3, and //y ex //x oc k~'^ are strongly irrelevant (in the 
renormalization group sense). Introducing the effective 
dimension D = d + 1, one finds h{r,t) ex k*^^"^^/^ and 
h{r,t) ex k'-'^+'^'/^. For the nonlinear couplings, we ob- 
tain A^ ex K^'^~^y^ and A|| ex k~^/^. Since D is clearly 
positive, the coupling A|| becomes irrelevant. The upper 
critical dimension dc for the theory is determined by A^ , 
via = 4— D which leads to dc = 3. The invariant dimen- 
sionless effective expansion parameter is T|| X^k^'^'^^^ 

as shown by the rescaling A_l -^ a'^/^Aj_, T|| —> a^T||. 
At the Gaussian level, this case corresponds to a critical 
line parametrized by Ty . The mean-field values for the 
roughness exponents are simple: In real space, one has 
X|| = (3 — d) /4 and x± = (3 — d) /2 while the momen- 
tum space exponents are given by X|| = (2 — d)/2 and 
X± = (4 — d) /2. While the potentially negative value of 
X|| might appear startling, it is simply a consequence of 
forcing Eq. (^J into the form of Eq. ^^. 

Case ii: Here, tj_ remains finite while tj| vanishes. The 
Gaussian part of the functional is dominated by hd^M 

and T^hV\h. Again, even at the tree level, parallel and 
transverse momenta scale differently: now, g|| ex k and 
|q^| ex K^+'^ with A = 1. Time scales as k~^. The 
strongly relevant perturbation is tii ex k^. One may 
still write h{r,t) ex k(^-4)/2 and h{r,t) ex k('D+4)/2 



but the appropriate effective dimension is now D = 
(rf— 1)(1 + A) + 1. The two nonUnearities switch roles so 
that A|[ ex k(s^^)/2 and Aj^ ex k^^/"^. In this case, A^, 
/Li_L, and /ix are irrelevant while Aii becomes marginal 
at the upper critical dimension dc = 7/2. The invari- 
ant dimensionless effective expansion parameter follows 
from the rescalings as /X|| t_|_ A||K^^''~^^/^. Again, 

it appears as if this case corresponds to a critical line 
parameterized by t^. However, we will see below that 
the order of the transition may well become first order 
once fluctuations are included. We therefore refrain from 
quoting mean-field roughness exponents here. 

Case in: Finally, we consider the multicritical point 
where both t±_ and t\\ vanish. Both momenta scale iden- 
tically, as q\\ ex Iqj^l ex k, so that A = at the tree level. 
We choose (3 so that /x^ scales to 1. One obtains ^t ex 
K~^, r|| ex r^ ex k^, and h(r,t) ex k^-'^"'*^/^, h{r,t) ex 
^(■D+4)/2^ with D ^ d. Both nonlinear couplings. Ay 
and X±, have the same upper critical dimension dc = 6. 
The effective expansion parameters are w = fj,x/^/JJ^, 

U|| = M^'/''A,|At('^-6)/2, ^^^ ^^ = ^-3/«AiK('^-6)/2^ j^ 

this case, the anisotropy exponent A vanishes at the tree 
level so that, to this approximation, all four roughness 
exponents are equal, given by (4 — d) /2. We will see, 
however, that this changes already in first order of per- 
turbation theory. 

In the following, we compute the scaling properties of 
correlation and response functions for the physically most 
interesting case (i) in a one-loop approximation. Our 
findings for cases (ii) and (iii) are reviewed only briefly, 
leaving the full technical analysis toi^. 



B. The one- loop approximation. 

We use dimensional regularization combined with min- 
imal subtraction22iSS. The basic building blocks of the 
perturbative analysis are the one-particle irreducible ver- 
tex functions F^^ jy({q, lj}) with N (TV) h (h-) ampu- 
tated legs. The notation {q, oj} is short-hand for the 
full momentum- and frequency-dependence of these func- 
tions. Focusing on the ultraviolet singularities, only those 
F^ jY with positive engineering dimension need to be con- 
sidered. Taking into account the symmetries and the 
momentum-dependence carried by the derivatives on the 
external legs, the set of naively divergent vertex func- 
tions is reduced to Fi^i and Fi_2 which are computed to 
one-loop order. Some technical details are relegated to 
the Appendix. 



Case i: t± 



and r« > 



This is the simplest non-trivial case. Only one param- 
eter, r_L, needs to be tuned to access criticality. Since A|| 
is irrelevant, it may be set to zero. Neglecting all other 



irrelevant terms as well, the functional simplifies to 

J[h,h] - -f fd''xdt{h[j-'dt + iVlf-Tiidl-T^Vl]h 
- h^ ~ X^h{diih)Wlh\ . (17) 

Thanks to the momentum dependence of the nonlin- 
ear term, X±h{d\\h)'S/']_h, all divergences in Fi^i and Fi^2 
are already logarithmic and appear as simple poles in 
e = dc — d. In a minimal subtraction scheme, we focus 
exclusively on these poles and their amplitudes to extract 
the renormalizations. Since the nonlinearity is cubic in 
the field, the expansion is organized in powers of A^^; i.e., 
the first correction to the tree level is always quadratic 
for Fi^i and cubic for Fi^2- The tilt invariance leads to 
a Ward identity connecting Fi^i and Fi^2, so that only 
the divergences in Fi^i need to be computed explicitly. 
Specifically, the tilt transformation h{r, t) — > h{r, t)+h-r, 
T± — > Tj_— ti||A_L shows that the parameter b renormalizes 
as the field h itself. Hence, the term X±h is renormalized 
by the same factor as tj_. 

Considering the perturbative contributions to 
Fi^i(q, w) further, we note that all of them carry 
external momenta, indicating that the terms hdfh and 
7/1^ do not acquire any corrections. Moreover, particle 
conservation in conjunction with invariance under 
parallel inversion and transverse rotations prevents the 
emergence of corrections to the T±hW\h term, at any 
order in perturbation theory. Hence, we should expect 
only two nontrivial renormalizations in Fi.i, namely 
those for the field h and for the parameter r|| . Leaving 
the detailed calculations of the renormalized quantities 
to the Appendix, we seek an infrared stable fixed point 
for the effective dimensionless coupling u, defined as 



A: 



-3/4 - 
-Til K 



-/^A 



(18) 



where A^ is a simple geometric factor, defined in the 
Appendix. A careful analysis of the flow equations for 
the renormalized parameters reveals the presence of a 
single infrared stable flxed point, at 



= ±4./-(l 



0{e)) , 



(19) 



where the sign is given by the sign of the initial coupling 
constant A^ . To obtain the scaling properties of correla- 
tion and response functions, we exploit the fact that the 
bare theory is independent of the external momentum 
scale K. The resulting partial differential equation (the 
renormalization group equation) and its solution is dis- 
cussed in the Appendix. It predicts, speciflcally, the scal- 
ing form of the height-height correlation function, C(r, t). 
Including the critical parameter, t±, in its list of argu- 
ments, we find: 

C(r,t;ri) = l'^+^-^+'JC{l'+^xi\,lr^,lH;l-^^''^T^) (20) 

where / is an arbitrary flow parameter. This behavior 
of C(r, t) is clearly a natural anisotropic gcneralizationi^ 



of the usual scaling form of critical dynamics. More- 
over, it is clearly consistent with the anticipated scaling 
behavior, Eq. Ullfl . especially if we set t± = 0. The ex- 
ponents v^, 7], and z have their usual meanings: v^ con- 
trols the scaling of the strongly relevant coupling tj_, 77 is 
the anomalous dimension of the field and controls criti- 
cal correlation functions, and z is the dynamic exponent, 
relating spatial and temporal fluctuations at criticality. 
A is the strong anisotropy exponent introduced in Sec- 
tion III. Exploiting the symmetries of the theory fully 
(see Appendix), we find that v±, rj, and z are related by 
scaling laws which are exact (at least within perturbation 
theory) : 



3. Case Hi: 



and T± 







V 



and 



l/i^i 



V 



Assuming that u* remains nonzero at higher orders of 
perturbation theory, its flow equation gives us another 
exact scaling law, relating A and 77: 

A + ?; = 2 - d/3 

As a consequence, only a single exponent, e.g., i], has to 
be computed order by order in perturbation theory. Our 
one-loop calculation results in 



r7 = -2e/3 + 0(e2). 



(21) 



Now, all others follow from exponent identities which are 
exact, at least as long as there are no non-perturbative 
corrections. 



2. Case ii: ry — > and t± > 

This is the second non-trivial case. Neglecting irrele- 
vant terms, the dynamic functional simplifies to 



J[h, h] =7 / d'^x dt [h[j-'dt + Miiaf - r||9f - r^Vi 



lhid\lh)d^h} 



(22) 



and the upper critical dimension is now dc = 7/2. 

Following the same approach as in the previous case, 
we first seek an infrared stable fixed point. The existence 
of such a fixed point guarantees that the model exhibits a 
scale-invariant regime where roughness exponents can be 
defined. As before, we define a suitable effective coupling 
M ex /in T^ A||k'-^'*^^-'/^ and analyze its flow. In 

this case, however, we find no stable fixed points corre- 
sponding to physically meaningful (i.e., real) values of the 
coupling M, at least to this order in perturbation theoryi^. 
The lack of such fixed points often indicates a first order 
transition, but only a more detailed analysis of the under- 
lying mean-field theory or a careful computational study 
will resolve this issue. 



Finally, we briefly review our results for the multicrit- 
ical point where both critical parameters vanish simul- 
taneously. This case was previously studied by Marsili 
et ali^, in a momentum shell decimation scheme. In 
this procedure, a hard momentum cutoff prevents the 
emergence of ultraviolet divergences in the momentum 
integrals. However, in theories with strong anisotropy 
and nonlinearities carrying multiple derivatives, the per- 
turbative corrections depend on how the cutoff is im- 
plemented, requiring extreme care. This may explain 
why we were unable to reproduce the earlier— results. In 
contrast, our findings are easy to check, since our field- 
theoretic approach does not suffer from these complica- 
tions. 

The full functional, defined by Eqs. Hl()|) and lO, now 
comes into play. To ensure the stability of the critical 
theory at the tree level, we demand /iy gy -I- 2/ixgjj <i± + 
(qi)^ > 0. This limits the physical range of fi\i and Hx to 
/i|| > and fix > ^\/JM\- To complicate matters further, 
both nonlinear couplings. Ay and Aj^, are marginal at 
the upper critical dimension dc ^ 6. While the detailed 
calculations become more involved, the technical analysis 
remains straightforward^^. In particular, thanks to the 
Ward identity, all renormalizations can still be obtained 
from the two-point function, Fi^i. 

The key results are as follows. Again, we define appro- 
priate effective couplings. One of these. 



w : 



M> 



appears in the propagator and generates w-dependent 
algebraic coefficients in the perturbation expansion. The 
other two, 



Mil 
U I 






control the nonlinar terms and are treated order by or- 
der in perturbation theory. C^ just absorbs some com- 
mon geometric constants. A careful analysis of the fiow 
equations for these three couplings showsi^ that there is 
only a single, physically meaningful infrared stable fixed 
point, given by 



(23) 



This contradicts earlier results^ where a fixed point with 
u*^^ was supposedly found. 

Exploiting the symmetries of our theory fully and as- 
suming that u*^ = remains valid to all orders in per- 
turbation theorjii^, we can determine the critical expo- 
nents associated with this model. Remarkably, we find 




that only a single exponent must be computed explicitly 
within the e-expansion, e.g., A. The scaling behavior of 
the height-height correlation function obeys the general 
form given in Eq. (|2U|I . except that both critical param- 
eters now appear, each with its own scaling exponent, 
vj^ and i/|| , respectively. Again, this confirms the antic- 
ipated scaling, Eq. Hll|l . The exponents, however, take 
different values here, demonstrating that cases (i) and 
(iii) fall into distinct universality classes. Our one- loop 
calculation yields 



(24) 



At the multicritical point, i.e., case (iii) with t±_ 
0, Til — > 0, we found 77 = and a nontrivial A. Hence, 



^ 23 + 6V15 ^, 2n 



The remaining exponents rj, z and v^ do not acquire any 
corrections beyond the tree results, so that 

'7 = 0, 2 = 4, ^-L=2 

to all orders in perturbation theory. The exponent vn is 
related to A through an exact scaling relation, namely. 



'W 



d-2+3A 



(25) 



These results are sufficient to evaluate the associated 
roughness exponents. 



V. RESULTS FOR THE ROUGHNESS 
EXPONENTS 

To begin with, we recall that the scaling forms of the 
two-point correlations for cases (i) and (iii) are indeed 
consistent with Eq. (|11|) . We may therefore immediately 
express the roughness exponents, Eqs. H13|) and H15|l . in 
terms of the exponent A and 77 for the two universality 
classes. 

For case (i), characterized by r^ — *■ at positive 
Til, we found only one independent exponent, namely 
T] = -2e/3 + 0(£2) with A = 2 - d/2 - 77. Writing 
all four roughness exponents in terms of 77, we arrive at 
expressions which are exact to all orders in e = d ^ 3: 



X-L = 1 - cf/3 x\\ 

X± = (4 - rf - 77)/2 XII 



l-rf/3 
3 - d/3 - 77 

4-77 d 

3 - d/3 - 77 ^ 2 



The mean-field values are easily recovered by setting e = 
0. The physically most interesting case corresponds to a 
surface grown on a two-dimensional substrate, i.e., d = 2 
and £ = 1. For this situation, one obtains x_l = 1/3 
and XII = 1/9 + O(e^) while X|| = -2/9 + 0{e'^) and 
X±_ = 4/3 + 0(£^). Remarkably, the exponent xj_ = 1/3 
is actually exact, at least to all orders in perturbation 
theory. 



X± = ^[4-(d + A)] XII 



Xx = -(4-d) X|| 



4 - (d + A) 
2(1 + A) 
2 d 

1 + A ~ 2 ■ 



All of them are negative near the upper critical dimension 
dc = 6. In order to access the physical {d = 2) situation, 
one has to set £ = 4 here which gives a huge anisotropy 
exponent, A ~ 16.814, if one naively uses the one-loop 
result. While roughness exponents can in principle be 
calculated, we should not place much confidence in their 
numerical values. Nevertheless, we still obtain testable 
predictions for this case, namely, the general scaling form 
of the height-height correlations as well as the scaling 
laws relating different exponents. 

Needless to say, in the absence of an infrared stable 
fixed point, the renormalization group gives us no infor- 
mation about possible roughness exponents for case (ii) 
{t± > 0, r|| -^ 0). In fact, if the scenario of a first or- 
der transition line applies, the whole concept of scaling 
exponents would be misplaced. 



VI. CONCLUSIONS 

To summarize, we have analyzed the long-time, large 
distance scaling properties of a surface growing un- 
der ideal MBE-type conditions, subject to an inci- 
dent particle beam tilted away from the surface nor- 
mal. This selects a particular ("parallel") direction in 
the substrate plane, so that the resulting growth equa- 
tions for the (single-valued) height field arc spatially 
anisotropic. In particular, an effective surface tension be- 
comes anisotropic, contributing the terms T||9j?/i-f rj^V^/i 
to the right hand side of the Langevin equation. If both 
of these parameters, th and tj^, are positive, the surface 
is described by the Edwards- Wilkinson model. However, 
depending on experimental control parameters such as 
temperature, incident flux rate and angle, or particle size, 
cither t\\ or t_l or both can vanish, generating signifi- 
cantly different surface properties. Each of these three 
possibilities leads to a distinct field theory with different 
upper critical dimension. We find dc = 3 if rj^ ^ at fi- 
nite T|| , dc ~ 7/2 if r|| — + at finite r^, and dc = 6 if both 
vanish. Only the first and the third case lead to scale- 
invariant behavior; the second one may in fact trigger a 
first-order phase transition whose properties lie outside 
the scope of our RG techniques. Focusing on the first 
and third case, we find two distinct, novel surface uni- 
versality classes. For both, we compute the scaling be- 
havior of the two-point height-height correlation function 
and carefully extract four different roughness exponents. 
Two of these characterize the height fluctuations of the 
surface in real space, scanned either along the parallel or 
the transverse directions; the remaining two character- 
ize scattering data with momentum transfer either along 



the parallel or the transverse directions. When analyzing 
experimental data, care must be taken in identifying the 
correct exponent. 

Clearly, the third case requires the careful tuning of 
two parameters, T|| and t±. If we can substantiate the 
presence of a first-order line for ry — *■ at finite t±, the 
point r|| ~ T± ~ would in fact be a critical endpoint, 
since it separates a line of second order transitions from 
a line of first order ones. In order to access it in a typ- 
ical experiment, at least two control parameters have to 
be tuned very carefully, making it difficult to observe. 
For this reason, either the second or the first order lines 
should be more accessible experimentally. From the RG 
perspective, even if ry is set to zero initially, it will be gen- 
erated under RG transformations, resulting in a non-zero 
value; this is not the case for the t_lV5^/i contribution. In 
that sense, we believe that the most physically relevant 
theory (apart from Edwards- Wilkinson behavior) is the 
one with r|| > and t± = 0. Amongst our key results 
for this model are the roughness exponents for real-space 
surface scans. For the physically most interesting case of 
a two-dimensional surface, we find X-L = 1/3 if the fluc- 
tuations are measured along the transverse direction, and 
X|| = 1/9 -|- O(e^) for scans along the parallel direction. 
While the value for X|| niay be modified by higher-order 
contributions in perturbation theory, the result for x± is 
exact. 

Finally, we turn to possible experimental evidence for 
these exponents. Thin films of vapor-deposited gold on 
smooth glass surfaces^Si^ were previously proposed^ as 
possible realizations of our theory. The deposit surfaces 
were imaged by scanning tunneling microscopy (STM), 
and their mean-square width was measured as a func- 
tion of STM scan length, resulting in a roughness ex- 
ponent of 0.35^2i2i. The temperature is sufficiently low 
(T = 298 K) so that desorption is negligeable. Even 
at small angles of incidence (between 2° and 25°), the 
growth is anisotropic. The STM images^ show some evi- 
dence for striped pattern formation, if larger sample areas 
(~ 3 X 10'^ nm on each side) are imaged. It is encouraging 
that a growth exponent close to 1/3 is observed since this 
matches our prediction for the transverse direction; un- 
fortunately, there seems to be no evidence for the much 
smaller exponent (1/9) which should control the parallel 
direction. A final resolution of these issues has to await 
a more detailed analysis of the experimental data. 



APPENDIX 

In the following, some technicalities associated with 
case (i) are presented. We will switch freely between the 
(q, t) and the (q, w) representations, depending on which 
one is more convenient. We first collect the elements of 
perturbation theory and then discuss the one-loop cor- 
rections. Neglecting all nonlinearities in Eq. H17|) allows 
us to identify the bare propagator, Go{q, t) , and the bare 
correlator. Co(q, t) , via 



Go(q,i)<5(q-q') 
Co(q,t)5(q-q') 



(/i(q,f)/i(-q',0))o 
e{t) exp(-r(q)7t) , 

(/i(q,i)/i(-q',0))o 
r(q)-iGo(q,|i|), (28) 



where 



r(q)=^||q|| + T||(7f|+T 



-qi 



Here, (•)o denotes a functional average with Gaussian 
(Aj_ = 0) weight. The Heaviside function 9{t) is defined 
with 6{0) = 0. Turning to the interaction terms, it is 
convenient to rewrite them in a symmetrized form. In 
Fourier space, the expression for the three-point vertex 
reads 

i^(qi,q2,q3) = «7A± qi_L-('72||q3±+q2_L'73|i)-'?i||q2_L-q3± 

(29) 
with qi + q2 + q3 = 0. The sign convention is such that 
all momenta attached to a vertex are incoming. 

We are using dimensional regularization, so that ultra- 
violet divergences appear as simple poles ui e = dc — d, 
with dc = 3. In a minimal subtraction scheme, we fo- 
cus exclusively on these poles and extract the renormal- 
ization constants from their amplitudes. Thanks to the 
symmetries of the theory, we anticipate only two nontriv- 
ial renormalizations, namely those for the field h and for 
the parameter rii , both of which can be obtained from 
Fi 1. To one loop order, the explicit expression for the 
singular part of the two-point vertex function is given by: 

ri,i(q,w)poic = *w + 7[T||qff + (qi)2+Tiqy (30) 

+ 7g[2r||gf-(qi)^]+0(.^), 

where the effective expansion parameter u is given by 
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A.r,r^/*«-^/2A, 



(31) 



and Ag summarizes a geometric factor which appears in 
all Feynman diagrams: 



''~\V^t('-^]t^' + ' 



^'~ (27r)'''"^ V 2 y ^ V 2 
Sd is the surface area of the d-dimensional unit sphere. 
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Keeping in mind that there are only two nontrivial 
renormahzations to all orders, we introduce renormalized 
quantities via 



A 



Tl = Z Tl 



'^11 — Z ZrT\l 



A 



Z-3/2A, 



The renormalized vertex function Fi i is defined by de- 
manding that 

^i,i{(Ii^i1,T-L,ti\,u,k) = ri,i(q,w,7,fj_,f||,A±) (32) 
be pole-free. One finds 

Z=l + — + 0{u*), 

oS 
2 



Zr = l 



4e 



0{u^ 



(33) 



The corresponding Wilson functions are defined as the 
logarithmic derivatives of the associated Z-factors, at 
constant bare quantities, i.e.. 



Asymptotic scaling results from this partial diff'crcntial 
equation at an infrared stable fixed point that is a solu- 
tion of (3{u*) = with f3'{u*) > 0. The one-loop approx- 
imation Eq. H35|l leads to the stable fixed point 



±4,/-(l + 0(£)) , (37) 

where the sign is given by the sign of the initial coupling 
constant A^. Thus, u* is non-zero in the e-cxpansion, 
and making the reasonable assumption that this remains 
true to all orders in perturbaton theory, we find from 
Eq. H35|) the exact relation 



C+Cr^Y 



(38) 



where we have denoted the Wilson (^-functions, evaluated 
at the fixed point, by a superscript *. Eq. H36|l can be 
solved easily at the fixed point u*, using the method of 
characteristics, combined with dimensional analysis and 
the rescaling invariances. If wc suppress unnecessary ar- 
guments, the solution can be written in the form 



C = tid^\nZ\ 



bare 



+ 0(u4) , 



Cr = Kdi^\nZr\hmc ^ ^ + Oiu'^) 



(34) 



The fiow of the dimensionless effective coupling constant 
u under renormalization is controlled by the Gell-Mann- 
Low function, 



r £ 3 



= "L-2 + 32" 



0{u*)]. 



(35) 



The renormalization group equation (RGE) for the 
Green functions, i.e., the connected correlation functions 
with N {N) external h {h) legs, simply states that the 
bare theory is independent of the external momentum 
scale k: 



dn 



= K-j--G^jj{{r,t};T±,Tii,'j,X±_) 



= ^— Z(^-^)/2G^,^({r, t}; r^, T||; u, 7, n) . 

In our case, this equation takes the form 
d ^ d ^ d ^ d 

OK OU 0"f OT± 



where I is an arbitrary fiow parameter. The two criti- 
cal exponents A and 77 can be expressed in terms of the 
Wilson functions as 



A = l 



c - c* 



v^C 



(39) 



while scaling laws give us the remaining exponents v, 
and the overall scaling exponent 6j^ f^: 



4-77, 

N 






(40) 



^N,N 



-(d + A-4 + 77) + -(d + A + 4-77) 



Finally, Eq. (|38|l provides another exact scaling law re- 
lating A and 77, provided u* ^ holds to all orders: 

A + 77 = 2 - d/3 



As a consequence, only a single exponent, e.g., 77, has to 
be computed order by order in perturbation theory, and 
all others follow from exponent identities. 
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